Raman scattering near a quantum critical point 
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Electronic Raman scattering experiments in a wide variety of materials (ranging from mixed-valence 
materials to Kondo insulators to high-temperature superconductors) show anomalous behavior in 
the -Bi 3 channel when the system is on the insulating side of the metal-insulator transition. Here we 
provide an exact solution for Bi g Raman scattering in the Falicov-Kimball model, show how these 
theoretical results are universal near the metal-insulator transition and show how they produce the 
two main features seen in experiments near a quantum critical point: (i) the rapid appearance of 
low-energy spectral weight as T is increased from and (ii) the existence of an isosbestic point 
(where the Raman response is independent of T at a characteristic frequency). 



Raman scattering has long been applied to study met- 
als, insulators, semiconductors, and superconductors. 
Via light's coupling to the electron's charge, inelastic 
light scattering can reveal electron dynamics over a wide 
range of energy scales and temperatures. More recently, 
Raman scattering has played a strong role in elucidating 
the nature of electron dynamics in the high temperature 
superconductors over the entire region of their phase di- 
agram. These studies have been particularly useful in 
shedding light on superconducting and pscudogap ener- 
gies and on magnon scales and dispersions. Systems as 
disparate as mixed-valence compounds (such as SmBg 
JjJ), Kondo- insulators (such as FeSi §), and the under- 
doped cuprate high temperature superconductors |j|-||), 
show temperature-dependent B\ g Raman spectra that 
are both remarkably similar and quite anomalous, sug- 
gesting a common mechanism governing transport. As 
these materials are cooled, a pile up of spectral weight 
appears for moderate photon energy losses with a simul- 
taneous reduction of the low frequency spectral weight. 
This spectral weight transfer is slow at high temperatures 
and then rapidly increases as temperature is lowered to- 
wards a putative quantum critical point (corresponding 
to a metal-insulator transition). In addition, the spec- 
tral range is divided into two regions: one where the 
Raman response decreases as T is lowered and one where 
the response increases. These regions are separated by 
a so-called isosbestic point, which is defined to be the 
characteristic frequency where the Raman response is in- 
dependent of temperature. 

While the existence of anomalous features in Raman 
scattering near a quantum critical point has been seen 
for some time, there is no theoretical understanding that 
connects the metallic and insulating states. In 1991, 
Shraiman and Shastry || outlined a procedure to con- 
struct a theory that can interpolate from a metal to an 
insulator. However, due to the ever increasing Hilbert 
space needed to describe a metal from the insulating side, 
or the lack of a clear picture of quasiparticles from the 
metallic side, quantitative calculations were not feasible. 
One is then left with approximate methods (such as per- 



turbation theory for the Hubbard model) or more phe- 
nomenological approaches to construct a Raman theory 
which suffer the limitations of being unable to reach dif- 
ferent phases and of lacking a microscopic basis. 

Here we provide the first exact theoretical descrip- 
tion of Raman scattering near a quantum metal-insulator 
transition that contains all of the anomalous behavior 
seen in experiments on these strongly correlated materi- 
als. We choose the spinless Falicov-Kimball model |7j as 
our canonical model for Raman scattering. It contains 
two types of electrons: itinerant band electrons and lo- 
calized (d or f) electrons. The band electrons can hop be- 
tween nearest neighbors [with hopping integral t*/(2y/d) 
on a d-dimensional cubic lattice] , and they interact via a 
screened Coulomb interaction with the localized electrons 
(that is described by an interaction strength U between 
electrons that are located at the same lattice site). We 
measure all energies in units of t* . The Hamiltonian is 
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+uj2 d UiWi, (i) 



where d\ (di) is the spinless conduction electron creation 
(annihilation) operator at lattice site i and Wi = or 1 
is a classical variable corresponding to the localized /- 
electron number at site i. We will adjust both Ef and 
\i so that the average filling of the (i-electrons is 1/2 and 
the average filling of the /-electrons is 1/2 (p, = U/2 and 
E f = 0). 

The Raman response is found from the frequency- 
dependent density-density correlation function that is de- 
picted in Figure 1. The Raman scattering process in- 
volves a two-photon-electron-hole vertex function that is 
called the Raman scattering amplitude 7(k). In addition 
to the effects of the Raman scattering amplitude, the 
density-density correlation function is, in general, also 
renormalized by the irreducible dynamical charge vertex, 
which is denoted by T in Figure 1. Evaluating the di- 
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FIG. 1. Dyson equation for the nonresonant Raman re- 
sponse function. Solid lines denote electron propagators and 
wavy lines denote photon propagators. The shading denotes 
the fully renormalized susceptibility and the symbol F is the 
irreducible frequency-dependent charge vertex. 
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FIG. 2. Bi g Raman response as a function of U at T = 0.5. 



agrams in the standard fashion produces the following 
result: 



k J ° 



(2) 



Tr(e-^p k (0)) 



with the uniform (q = 0) Raman density operator 

Pk = 7(k)4d k , d k = i^ e - R ^ k d i5 (3) 



Z = Tr(e~P H ) , the partition function, and ivi = 2iirlT 
the bosonic Matsubara frequency (the r-dependence of 
the operators is with respect to the full Hamiltonian). 
The Raman scattering amplitude 7(k) is a complicated 
function of the incoming and outgoing photon polariza- 
tions, of the photon energies, and the polarizability of the 
medium. In nonresonant Raman scattering one neglects 
the frequency dependence of the Raman scattering ampli- 
tude, and characterizes the Raman response in terms of 
the different spatial symmetries of the remaining function 
7(k). One can expand this function in a Fourier series 
and examine the contributions of the lowest components 
of the series, and compare them to experiment (the low- 
est order B\ g contribution is 7(k) = 1) J coskj, 
with d — > oo the spatial dimension). More sophisticated 
approaches would calculate the Raman scattering ampli- 
tude from "first-principles" and would include any pos- 
sible resonant Raman scattering effects. We leave those 
pursuits to future work. 

The Falicov-Kimball model can be solved exactly in the 
infinite-dimensional limit by using dynamical mean-field 
theory (see Ref. || for details). The dynamical charge 
vertex is local in infinite dimensions which implies that 
correlation functions that have the same symmetry as 
the lattice are renormalized due to this charge vertex, 
but correlation functions that are orthogonal to the lat- 
tice, have no vertex corrections, and so they are repre- 
sented by their bare bubble diagrams M. Although the 



fully symmetric Ai g Raman response can be calculated 
directly, we concentrate here on the B\ g response, which 
is described by the bare bubble diagram in Figure 1. In 
this straightforward symmetry analysis also shows 

that resonant Raman scattering effects vanish for nearest 
neighbor hopping on a hypercubic lattice in large dimen- 
sions, so the total Raman response is represented by these 
nonresonant results. 

The Falicov-Kimball model has a ground state that is 
not a Fermi liquid because the lifetime of a quasiparticle 
is finite at the Fermi energy. As U increases, the system 
first enters a pseudogap phase, where spectral weight is 
depleted near the chemical potential, and then under- 
goes a metal-insulator transition. The interacting density 
of states (DOS) is, however, temperature-independent 
for fixed U and fixed electron fillings. For half-filling, 
U < 0.65 corresponds to a weakly-correlated metal, while 
a pseudogap phase appears for 0.65 < U < 1.5 moving 
through a quantum critical point at U = 1.5 to the in- 
sulator phase U > 1.5 (we neglect all possible charge- 
density- wave phases here). 

The Raman response is calculated in a similar fashion 
to the dynamical charge susceptibility the Dyson 

equation of Figure 1 is formally analytically continued 
to the real axis, where the expression for the Raman re- 
sponse is found to depend solely on complex integrals of 
the interacting single-particle DOS. In Figure |^ we plot 
the Raman response at a fixed temperature T — 0.5 for 
different values of U. For small values of U, a small scat- 
tering intensity is observed due to the weak interaction 
among "quasiparticles" providing a small region of phase 
space allowable for pair scattering. The peak of the re- 
sponse reflects the dominant energy scale for scattering, 
as is well known in metals 12 and the high-energy tail 
is the cutoff determined by the finite energy band. As 
U increases, the low-frequency response is depleted as 
spectral weight gets shifted into a large charge transfer 
peak at a frequency ~ U . The charge transfer peak be- 
gins to appear for values of U for which the DOS is still 
finite at the Fermi level and becomes large in this pseu- 
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FIG. 3. Big (a) Theoretical Raman response as a function 
of temperature for U — 2 (which lies just on the insulating 
side of the metal-insulator transition) and (b) experimental 
Raman response for FeSi at moderate to low tempera- 
tures (where the isosbestic point and low temperature spec- 
tral weight depletion is seen). The experimental graphs are 
labeled by the temperature (in K) where the data was col- 
lected. 



dogap phase before growing even larger in the insulating 
phase. Notice how low-frequency spectral weight remains 
even as one is well on the insulating side of the quantum 
critical point. It is these spectral features that are char- 
acteristically seen in the experiments and which can only 
be seen in a theory that approaches the quantum critical 
point. 

In Figure ||(a) we plot the temperature dependence 
on the insulating side of the metal-insulator transition. 
For small values of U the spectra are only slightly de- 
pendent on temperature due largely to the small changes 
in kinetic energy with T. As U increases into the pseu- 
dogap and insulating phases, nontrivial temperature de- 
pendences begin to appear. The total spectral weight 
increases dramatically with decreasing temperature as 
charge transfer processes become more sharply defined. 
At the same time, the low-frequency response depletes 
with lowering temperatures, vanishing at a temperature 
which is on the order of the T = insulating gap. This 
behavior is precisely what is seen in experiment B on 
FeSi at low temperatures ||(b) where both the isosbestic 
point and the low temperature spectral weight depletion 
can be seen. 



We attribute the presence of a low-frequency response 
in a system which is a strongly correlated insulator to the 
appearance of thermally activated transport channels. In 
the insulating phase at zero temperature, the only avail- 
able intermediate states created by the light must involve 
double site occupancy of a conduction and a localized 
electron, with an energy cost of U. This gives the large 
charge transfer peak at an energy U. As the temperature 
is increased, for half filling, double occupancy can occur 
and as a result light can scatter electrons to hop between 
adjacent unoccupied states either directly or via virtual 
double occupancies. The number of electrons which can 
scatter in this fashion increases with increasing tempera- 
ture, leading to an increase in the low-frequency spectral 
weight. The frequency range for this low- frequency Ra- 
man response is determined by the lower Hubbard band- 
width, which is typically much larger than the tempera- 
ture at which these features first appear. If one were to 
interpret the temperature at which the Raman spectral 
weight starts to deplete as the transition temperature 
T c and the range of frequency over which the weight is 
depleted as the gap A, then one would conclude that 
near the quantum critical point 2A/fesT c 3> 1. This is 
because the "T c " is determined by the gap in the single- 
particle density of states (which approaches zero at the 
quantum critical point), while the "A" is determined by 
the width of the lower Hubbard band (which remains fi- 
nite at the quantum critical point); hence the ratio can 
become very large near the quantum critical point. 

The spectral weight transfer from low frequencies to 
the charge transfer peak as a function of temperature can 
be quantified by separating the Raman response into two 
regions determined by the isosbestic point and plotting 
the total low-frequency spectral weight versus tempera- 
ture (not shown). Choosing U/2 as the location of the 
isosbestic point, we find that the reduction of spectral 
weight from high to low temperatures is over 50 percent 
even in the weak pseudogap phase, and decreases by well 
over three orders of magnitude as U increases into the 
insulating phase (U = 4). 

Lastly, we plot the inverse slope of the Raman response 
in Figure |I](a) as a function of temperature for different 
values of U and the experimental || results on optimal 
and underdoped copper oxides ^(b). Since the self energy 
is temperature independent, we might expect a constant 
Raman slope as a function of temperature, as is the case 
with disordered noninteracting electrons. However, this 
is not the case due to the formation of a thermally gener- 
ated band for scattering. For small values of U, the tem- 
perature dependence of the Raman inverse slope is weak 
due to the temperature independence of the self energy. 
However, as the single-particle bands begin to separate, 
the relevance of thermally generated double occupancies 
becomes more pronounced and the inverse slope becomes 
temperature dependent at low temperatures. We see that 
as U increases the low temperature slope increases dra- 
matically due to the depletion of low-frequency spectral 
weight. In particular, even in the pseudogap phase the 
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FIG. 4. (a) Theoretical inverse slope of the Bi g Raman re- 
sponse and (b) experimental Q inverse slope for optimal and 
underdoped high-temperature superconductors (the curves 
are labeled by their superconducting transition temperature 
and whether they are under or overdoped; the 57K and 81. 7K 
data have been multiplied by constant factors to separate 
them) . 



Our theoretical results compare quite favorably to the 
experimental results seen in a wide range of different 
materials ranging from mixed- valence compounds jlj] , to 
Kondo insulators j^j to the underdoped high-temperature 
superconducting oxides P-^|. In particular, all of those 
experimental systems appear to be close to, but on the in- 
sulating side of the metal-insulator transition, and hence 
they illustrate the two characteristic behaviors seen in 
our theory: (i) there is a rapid rise in the low-frequency 
spectral weight at low temperatures (at the expense of 
the high-frequency spectral weight) and (ii) there is an 
isosbestic point. Our model always produces an isotropic 
gap, so we are unable to illustrate some of the behavior 
seen in the copper-oxides where only the B\ g response 
is anomalous, and the A± g and Bi g responses are metal- 
lic rather than insulating. But our results do indicate 
a "universality" and model independence of the Raman 
response on the insulating side of, but in close proximity 
to, a quantum critical point. We believe this is the reason 
why so many different materials show the same generic 
behavior in their electronic Raman scattering. 
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inverse slope rises with decreasing temperature indicat- 
ing the proximity to the quantum-critical metal-insulator 
transition. As U increases into the insulating phase, the 
temperature dependence of the Raman inverse slope is 
indicative of the formation of gapped excitations. Note 
how similar the experimental results are to the theoretical 
predictions (the lowest temperature data for the under- 
doped case (top curve) has large error bars because the 
signal is so small; in all cases as the system becomes un- 
derdoped the Big response behaves more insulator-like). 

In more complicated correlated models of the metal- 
insulator transition the single-particle density of states 
will have a Fermi liquid peak at low frequencies which 
may add new features to the Raman response, but on 
the insulating side of the transition, where most of the 
anomalous behavior is seen, the single-particle density 
of states should be very similar to that of the Falicov- 
Kimball model (except for some additional weak tem- 
perature dependence of the interacting DOS), which is 
why these results are generically expected to be model- 
independent. 
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